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PREFACE 



The research reported here was done at the U# S. Naval Postgraduate 
School, Monterey, California during the period February through April, 
1954* The work was undertaken to determine experimentally the maximum 
stress and the torsional rigidity of a round shaft with diametrically 
opposite parallel flat sides. This information was obtained for several 
configurations by direct experiment and for one configuration by an 
electrical analogy# 

The author is greatly indebted to Prof# R# E# Newton for many 
helpful suggestions and much constructive criticism and assistance in 
the preparation of this manuscript# In addition the author wishes to 
thank Mr# Norman talker for the extreme care and patience exercised 
in the preparation of the test specimens# It is further desired to 
acknowledge the help of F. N* Lauck, Ph3, USN, who did the photographic 
work which appears in the thesis# 
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SUMMARY 



This thesis describes and reports the results of an experimental 
determination of the maximum shearing stress and twisting moment of cir- 
cular shafts with diametrically opposite flat sides# Five AL-61 Aluminum 
alloy specimens were tested, one fully round shaft and the remaining four 
with progressively larger flats# Strains were measured with resistance 
wire strain gages. Angles of twist were determined by optical means# 

An electrical analog experiment was conducted to determine independ- 
ent results for one of the configurations. 

The experimental results obtained by both techniques compared very 
favorably with mathematical solutions obtained by previous investigators# 
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CHAPTER I 



AN INTRODUCTION TO THE PROBLEM 



1. Statement of the problem . 

This thesis is a report of an experimental determination of the 
torsional rigidity and maximum torsional stress occurring in a solid 
round bar having diametrically opposite flat sides. 

2 . State of stress for ji solid round shaft with diametrically opposite 
parallel flats . 

The state of stress can best be visualized by reference to the 
membrane analogy for the torsion problem Considering the de- 

flected surface of the membrane conclusions can be made as to the stress 
distribution. The shear stress in tiny direction is proportional to the 
slope of the surface perpendicular to that direction. It varies from 
a maximum tangent to the contours to zero perpendicular to the contours. 
Thus the maximum shear stress acts at points where the contour lines are 
closest to each other. With this in view reference to Fig. 1, which is 
a cross section of the shaft, leads ^ 

one to conclude that the maximum 
stress will be at the center of 

, 

the flat sides and zero at the 

center of the cross section and 

at the four comers. 



Fig. 1. Cross section of 
the shaft . 



3. Previous solutions to the problem. 




solution for maximum shearing 
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stress employing a variation of the teohnique given by Shepherd £ 5 J • 
Qkubo obtained an approximate solution in the form of elementary functions 
•which satisfied the differential equation but only approximately satisfies 
the boundary conditions. YV. J. Carter and J. B. Oliphint £ 4 J obtained 
solutions for torsional rigidity and maximum shear stress by numerical 
relaxation. For very small flats Carter and Oliphint made use of a 
hydrodynamic analogy to obtain the maximum stress. The above authors 
carried out the computations for several configurations, but no experi- 
mental results were presented. A search of the available literature by 
the writer disclosed no further work on the problem. 

4. Electrical Analogy 

The problem oonsidered is that of twisting a prismatical bar with 
couples applied at the end cross sections (St.-Venant torsion £'?£)• 

It is assumed that the bar is homogeneous, isotropic and obeys Hooke's 
law. -u. 




St.-Venant showed that the only stress components different from 



zero are and (Fig. 2 ) and that they are functions of x and 

y only. L. Prandtl introduced the stress function <p= ‘PU.y-I 

defined by 




x 



u) 



Fig. 3. Cross Section of the Shaft. 

The stress function satisfies the partial differential equation 






■^2 = —ZG O' in R 



r 



( 2 ) 



and the boundary condition 



0 — O on C. 



( 5 ) 



G is the shear modulus for the material and -0- is the angle of twist per 
unit length# R is the region of the cross-section of the shaft and C is 
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its boundary curve (see Fig. 3* ) The twisting moment is found by 




(j) sx 



( 4 ) 



R 

To transform the problem, assume the stress function has the form 



0 = 0 '- ^ w+y) 



(5) 



in the region R of Fig. 3 and is equal to zero on the boundary C. Put- 
ting this in Eq. (2 ) gives 

v a 0 / - O (6) 



•which is Laplace T s equation to be satisfied in R# The boundary condition 
now reduces to 

0 ' = ex 2 +y 2 ) on C. (7) 

Thus the torsional problem reduces to obtaining a solution of Laplace's 
equation over the cross section and satisfying the boundary condition 

Eq. (7). 

There are many techniques available for obtaining a solution of 
Laplace's equation. The author chose to use a General Electric Analog 
Field Plotter -with Teledeltos uniform conducting paper because of its 
simplicity, ease of use and the availability of components. The experi- 
mental assembly used is shown in Fig. 4. 
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Fig. 4. Analog Experimental Assembly. 
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<p'-K(V) 



( 6 ) 



whore K is a constant of proportionality and V is the potential# Then 
by Eq. (7) the potential varies parabolioally along the flat side and is 
constant along the arc# The equipotential lines are lines of constant 
(j)\ From the equipotential lines and K, the (j) contours can be obtained 
for the region. Knowing (^/enables one to obtain the (j) contours by 
means of Eq# ( 5 )• 

Eqs# (l) now become, by substitution of Eqs. (5) and (8), 



■>»=*$ - eey. , + <B-&x 



where and are the components of the potential gradient as 

dX 0 q 

obtained from the equipotential lines. Eq. ( 4 ) now becomes 



Mt=2KJ y V&X&y 




(x*-+y x ) ix Sy 



(4a) 



R R 

The first integral may be evaluated using the equipotential plot and the 
second analytically. Thus the complete state of stress and the torsional 
moment may be obtained from the electrical analogy. 
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CHAPTER II 

MATERIAL AND TESTING TECHNIQUE EMPLOYED 



5 • Specimen material and dimensions # 

Aluminum alloy A*L - 6l was chosen for specimen material because of 
its low shear modulus and ease of machining. A low shear modulus was 
desirable in order to produce relatively large strains and angles of 
twist of the specimen for a given torque, both of which tended to reduce 
the relative error of the results# A three and one-quarter inch diameter 
shaft was chosen# From Fig. 5 it can be seen that the stress pattern 
must change in going from the small chuck diameter to the relatively lar- 
ger section, this necessitated having a long specimen in order that the 
ends of the ten inch gage length be approximately three diameters from 
the section change# By applying Saint- Venant * s principle the stress 
condition in the gage length could be considered as uniform# 

6# Method of machining specimens . 

The specimens were machined in accordance with Fig. 5 and Table 1# 



TABLE 1 

SPECIMEN DIMENSIONS 



C < 


y ~ 


ult 


b 


0 ° 


1.627 in. 


3.255 in. 


0 


11.0° 


1.576 in. 


1.545 in. 


0.602 in. 


2 5.0° 


1.576 in. 


1.428 in. 


1.331 in. 


46.0° 


1.576 in. 


1.055 in. 


2.2 65 in. 


65.7° 


1.576 in. 


0.646 in. 


1.436 in. 
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The specimens were first machined round and then the flats milled. 

All the specimens were cheoked for axial looatian and parallelism 
of the flats and were found to be within one thousandth of an inch which 
was considered satisfactory. 

7 • Strain gage selection and application . 

The selection of the proper strain gage was an important considera- 
tion as it was to be used in the determination of the maximum stress on 
the flats. The strain varied from zero at the comers to a maximum at 
the center of the flats, this dictated the selection of a gage as small 
as practicable. The principal strain directions are at an angle of A'f 
to the axis of the bar for torsion so a single element gage was chosen. 
Due to the A’f orientation, the width of the gage element was equally as 
important a consideration as the length. Two gages were used, one on 
each flat oriented so as to be ninety degrees to eaoh other, this enabled 
the separation of unwanted bending strains introduced during testing. 
Baldwin SR - 4 strain gages were used. 

A-19 gages with a one- sixteenth inch by one-sixteenth inch gage 
element were initially chosen. The use of these gages was abandoned 
since consistent results could not be obtained when they were tested on 
the fully round specimen (°C=d). Table 2 gives the results obtained 
from an A-l? gage. 
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TABLE 2 



A-19 STRAIN GAGE RESULTS 



Date 


Total Drying Tine At 
Room Temperature 


Total Drying Time 
At 1 60-17 0°F 


Apparent 

Strains 




4 March 


72 hours 




0.012J yW 


in/in 

in-lb, 


6 March 


120 hours 




0.0154^ 


in/in 

in-lb 


12 March 


264 hours 


20 hours 


0.0150 yU 


in/in 

in-lb 


18 March 


408 hours 


62 hours 


0.01 54yw 


in/in 

in-lb 



Two separate A-19 gages were tried but the same phenomenon persisted* 
Fig* 6 is a typical curve obtained from an A-19 gage, the relaxation on 
unloading was prevalent on all runs* In view of this, A-7 gages were 
substituted and they gave consistent results* The A-7 gages have a 
one-fourth inch gage length and three thirty-seconds inch width. Two 
AX- 5 gages were installed on the fully round specimen for preliminary 
testing and calibration of the single element gages* All gages were in- 
stalled and dried using the technique recommended by the manufacturer* 
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6* Angle of twist measurement s 

The maximum angle of twist was less than ^ of 1 degree so an 



optical means was chosen to measure angles of twist* The schematic 
diagram for the optical arrangement used at each end of the 10 inch 
gage length is shown in Fig* 7* 



The essential elements are a hairline illuminated by a light source 
so arranged that the rays from the hairline are reflected by a mirror 
mounted on the specimen and focused by a galvanometer lens, on a frosted 
scale* The actual equipment used can be seen in Fig* 8* The hairline 
images from both the light sources were focused on one frosted scale* 
Thus differences in the positions of the images on the frosted scale was 
a measure of the angle of twist in the gage length* Since small angles 
were involved, the sine, tangent and angle in radians were considered 
equal* The formula for the angle of twist for a given load is. 




collimating lens 



Fig* 7* Schematic diagram for optical system used for angle of 



twist measurement* 



where A /“) is the change in distance between the hairline 
images and d is the mirror-to-scale distance. The angle of twist per 
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unit length for a given load is -Q- —~Jq 1, ^ ere 10 inohea is the gag© 
length. 

9 . Strain gage tests . 

The equipment as used is shown in Fig. 6. The speoimens were 
loaded in a standard Tiniua - Olson Torsion Tester on the 0 to 50,000 

Q 

in-lb range except for the thin specimen ( 0< =z 65/7) when the 0 to 
12,000 in-lb range was used# Two AX - 5 strain gages mounted on the 
fully round specimen at physical locations ninety degrees apart were used 
to determine the best method of putting specimens in the machine without 
introducing extraneous loadings. If the rotating head of the machine 
was at or near zero degrees and the chucks were tightened slowly and 
alternately the maohine applied essentially pure torsion. Two runs were 
made on each specimen, one with each gage in the up position. A Baldwin 
Southwark Type K Strain Indicator was used with a Baldwin S R - 4 12 

Channel Switching and Balancing Unit to measure all strains* The temp- 
erature compensating gages were mounted an a separate piece of the same 
material. The specimens were loaded and unloaded in equal increments. 
After eaoh increment of load had been applied the strains, hairline 
image positions on the frosted soale and load were reoorded. Strains 
were read to the nearest lt~l/in$ load to the nearest 10 inoh-lb. 
and hairline image position to the nearest tenth of a millimeter* 

10* Analog Experiment 

The equipment used for the analog experiment is shown in Fig. 4. 

It is a General Eleotric Analog Field Plotter. The experiment was car- 
ried out only for ©<— -^5 , ~W ~ 30 /ft • As can be seen in Fig. 5# 
©< is the angle between the x axis and the exterior oomer of the 
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Fig. 8. Experimental Equipment Set Up 

15 




shaft, and w is ^ the distance between flats. Due to the double symmetry 
of the cross section the conducting sheet represented only ^ of the cross 
section as shown by the shaded area in Fig* The large dimensions 
were chosen to increase the accuracy of the results* 

The null method of determining potentials made it convenient to 
record experimental potentials P on an arbitrary scale. Accordingly, 
the potential at the center of the flat (w ,0) was taken a3 0 and that 
on the curved portion of the boundary was taken as 1* As is ^hown later, 
the potential V of Eq* 8 is related to the experimental potential P by 
the equation V-l+P. 

The parabolic potential distribution along the boundary X — IaT was 
approximated by a series of linear potential drops# Fig* 9 shows the 
desired parabolic distribution and how it was approximated by a series 
of linear drops# There was no net gain or loss of area under the curve* 
Small increments were taken at the center of the flat in order to have 
the smallest deviation from the desired distribution, and keep errors in 
potential gradient to insignificant magnitudes in this region. Table 5 
gives the potentials maintained at the chosen locations along the boundary. 
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TABLE 3 



POTENTIALS MAINTAINED ALONG BOUNDARY 



Y (in. ) 

0 

1 

s 

2 
4 

7 

12 

17 

22 

27 

30 



p 

0 

0 

0.0039 

0.0162 

0.0496 

0.1554 

0.3165 

0.5332 

0.8054 

1.000 



The experimental potentials measured along the boundary are plotted 
in Fig. 9. To establish the above potential points along the boundary with 
linear drops between points, 0.004 inch diameter Advance wire with lead 
wires soldered to it at the selected locations was clamped onto the 
boundary of the paper by the clamping arrangement shown in Fig. 4. The 
lead wires were connected to a slide wire potentiometer which was across 
the output voltage of the rectifier unit. The actual boundary potentials 
were established 0.15 inch in from the wire by varying the positions of 
the lead wires on the slide potentiometer. It was found that if the 
correct potential distribution was set along the boundary wire the act- 
ual potential distribution in the adjacent paper was significantly dif- 
ferent. The constant potential along the curved portion of the boundary 
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was maintained by use of silver paint backed by a copper strip, 
equipotentials were plotted using standard null technique. 



The 
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CHAPTER III 



REDUCTION OF DATA 



11. Strain ^a 3 e calibration * 

The AX-5 strain gages on the fully round specimen were used to 

calibrate the single element A-7 gages* This was necessary because the 

gage factor given by the manufacturer gives correct results for uniaxial 

stress in the direction of the gage length* The gage factor for this 

strain condition was found from 

A-7 Manufacturers gage factor 

c A-7 gage factor for this strain condition 

where £ is the strain in the direction of the A-7 gage as obtained 

from the AX- 5 gages and g is the apparent strain obtained by the A-7 

gage when the manufacturer’s gage factor was set on the indicator. 

The manufacturer’s gage factor was 1*95 and the apparent ga^e factor for 



this strain condition was found to be 1*97* It is noted that this re- 
sult implies a negative transverse sensitivity for the A-7 gage. This 
is in agreement with transverse sensitivities reported by Baumberger 
and Hines f 7 J • 



i 
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12 . Electrical Analogy 



The computation of the maximum stress at the center of the flat and 
the computations for torsional moment v;ere carried out as discussed in 
section 4. By Eqs. (Q) and (5) ; (p'-K(V) - ( ~X x i- ^ ) 

At (- WjO ) ; KV, and at (w, uw) ; ft V z - (5 ■& W* • 
Taking the applied potential difference as unity, — 1/ = / 

Then K— — ^ V z ~ 2. ^ \Z t ~ I . The potential at any 

point on the surface is, V — 1/ + P - ! + P where P is the experi- 
mental potential. Fig, 10 shows the equipotentials over a quarter of 
the cross section. At (air. 0 J ~TZ? — 0 and 

> f-Jr 



-r - dtf' , r~ n - ' - 

+6 ' ex dx 



-4- G-0--ur 



1P| 



— 6 - 0 - 



- (-Ur, O) 

is obtained from Fig. 11. Eq. (4a) now becomes 

r 



i 



(l+PlS)X&y. - e-6- J j(H x + y 3 -iS-X Sty (4a) 



Simpson’s one-third rule for approximate integration was used to obtain 
the volume under the potential surface for evaluating the first integral. 
The value of the last integral was obtained analytically. 



21 




Fig. 10. Equipotential Traces of Electrical Analog. 
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CHAPTER IV 



RESULTS AND DISCUSSION 

13* The results . 

Numerical results obtained from tests are given in Table 4. The 

o 

values for c< = 45 were obtained by the electrical analogy. 

TABLE 4 

VALUES OF MAXIMUM SHEAR STRESS 
AND TWISTING MOMENT 



CK 


TrnaX 




0° 


1.00 


1.58 


11° 


1.15 G -& -ur 


1.57 &■& r 4 " 


2 5° 


l.Jl G-0 -tw" 


1.44 


45° 


1.66 G-9-zo~ 


0.564 G-Or 4 * 


46° 


1.56 ( B ' O ' ~ ws ~ 


0.523 G'O-r 4 ' 


65.7° 


1.67 Q-Q-iat 


0.257 G£~r^ 



r is the radius of the shaft. 
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14. Comparison with previous solutions . 



Fig. 12 gives a comparison of the present results with those of 
previous investigators. Carter and Oliphint stated that they believed 
their results were accurate to within 2 or 5 per cent. 5kubo gave no 
figures for accuracy of his results. 
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15. Sources of error . 

The strain vs torque relationship obtained with the AX-5 strain 
gages an the fully round specimen was linear within reading accuracy. 
The extraneous bending stresses introduced by the torsion machine were 
found to be less than 1 % of the torsional stress. Results were cor- 
rected for bending stresses by averaging readings of gages on opposite 
flats. The orientation of the strain gages is believed to have been 
within 1° of the desired 45° orientation. Location of the mirrors on the 
10 inoh gage length was within - l/l6 inch. The strain gage elements 
covered a finite area on the flats of the shaft and the stress varies 
from zero at the exterior corners to a maximum at the center of the 
flat. Thus the strains indicated by the gages were not the maximum but 
seme sort of average over the gage element. This error is believed to 
be small. 

The manufacturer of the Teledeltos paper states that it has aniso- 
tropic variations of -5 to maximum deviations from a mean value of 
resistance. The rated linearity of the potentiometer used to establish 
the equi potentials was --g?o. No correction was made for the boundary 
potential being established 0.15 inch in the interior. 

The optical system used for measuring angles of twist gave reading 
accuracies of 1 % of the total angle of twist. The obtaining of greater 
acouraoy would have greatly complicated the optical system. 

The overall results are believed to be aocurate to within 4/£. 
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CHAPTER V 



CONCLUSION 

The experimental results checked closely with the previous 
mathematical solutions to the problem. Although the results obtained 
by the electrical analogy were reliable, the work involved was many 
times that of the direct experimentation for the same specimen 
configuration. 
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